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We have developed a general theory for a parametric pump 
consisting of a nonmagnetic system with two ferromagnetic 
leads whose magnetic moments orient at an angle with re- 
spect to each other. In this theory, the leads can maintain at 
different chemical potential. As a result, the current is driven 
due to both the external bias and the pumping potentials. 
When both and external bias are zero, our theory recovers 
the known theory. In particular, two cases are considered: 

(a) , in the adiabatic regime, we have derived the pumped 
current for arbitrary pumping amplitude and external bias. 

(b) . at finite frequency, the system is away from equilibrium, 
we have derived the pumped current up to quadratic order in 
pumping amplitude. Numerical results show interesting spin 
valve effects for pumped current. 
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I. INTRODUCTION 

Recently, tbfxp is considerable interest in the paramet- 
ric pumpingEI^ES. The parametric pump is facilitated by 
cyclic variations of pumping potentials inside the scat- 
tering syster&j and has been realized experimentally by 
Switkes et afl. On theoretical side, much progress has 
been made towards understanding of various features re- 
lated to the parametric-jDjuma. This includes quantiza- 
tion of pumped chargeaEln3'E3, the influence of discrete 
spatial symmetries and magnetic fielcB'LJ, the rectifica- 
tion of disalacement currentE£l, as well as the inelastic 
scatteringtil to the pumped current. The concept of op- 
timal pump has been raepposed with the lower bound for 
the dissipation derivedEJ. Within the formalism of time- 
dependent scattering matrix t hcoijy ■ , t hcJiaat current and 
shot noise in the pumping procesalZrEjoO has also been 
discussed. Recently, the original adiabatic pumping the- 
ory has been extended to account for the effect due to 
finite frequencylijtJ, Andijeev reflection in the presence 
of superconducting leacli3'EpL and strong electron inter- 
action in the Kondo regimecJ. This gives us more phys- 
ical insight of parametric pumping. For instance, the 
experimental observed anomaly of pumped current at 
(j) = and (j) — n can be explained using finite frequency 
theorylij as due to the quantum interference of different 
photon assisted processes. When superconducting lead 
is present, the interference between the direct refiection 
and multiple Andreev reflection gives rise to an enhance- 
ment of pump|ed current which is four times of that of 
normal systemll3. It will be interesting to further extend 
the parametric theory to the case where the ferromag- 
netic leads are present. With the theory extended, many 



new physics are foreseenEHl which may lead to pew oper- 
ational paradigms for future spintronic devicesE3. In this 
paper, we have developed a parametric pumping theory 
for a nonmagnetic system with two ferromagnetic leads 
whose magnetic moments orient at an angle 9 with re- 
spect to each other. Our theory is based on nonequilib- 
rium Green's function approach and focused on current 
perpendicular to plane geometry. Parametric pump gen- 
erates current at zero external bias. It would be inter- 
esting to see the interplay of the role played by pump- 
ing potential and external bias ifi-lihe leads are main- 
tain at different chemical potentialcj. Hence in our the- 
ory, the external bias is also included. In the adiabatic 
regime, the pumped current is proportional to pumping 
frequency. In this regime, we have derived the paramet- 
ric pumping theory for finite pumping amplitude. At 
the finite pumping frequency, the system is away from 
equilibrium, we have performed perturbation up to the 
second order in pumping amplitude and obtained the 
pumped current at finite frequencies. The newly de- 
veloped theory allows us to study the pumped current 
for a variety of parameters, such as the pumping ampli- 
tude, pumping frequency, phase difference between two 
pumping potentials, the angular dependence between the 
magnetization of two leads, as well as the external bias. 
We have applied our thejocy to a tunneling magnetore- 
sistance (TMR) junctionEiJ. Due to the reported room 
temperature operation of TMR, the fundmental principle 
and transport properiies of TMR devices has attracted 
increasingly attentior£3. Our numerical results show in- 
teresting spin valve effect for pumped current. The pa- 
per is organized as follows. In section II, we derive the 
general theory of a parametric pump in the presence of 
ferromagnetic leads. The numerical results and summary 
are presented in section III. 



II. GENERAL THEORY 

The system we examine consists of a nonmagnetic sys- 
tem connected by two ferromagnetic electrodes to the 
reservoir. The magnetic moment M of the left electrode 
is pointing to the z-direction, the electric current is flow- 
ing in the y-direction, while the moment of the right elec- 
trode is at an angle 9 to the z-axis in the x ~ z plane. 
The Hamiltonian of the system is of the following form 

H^Hl + Hr + Ho + Vp + Ht (1) 
where Hl and Hr describe the left and right electrodes 

Hl = Y^^ekL + aM) cl^^CkLa (2) 
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A. Pumping in the low frequency limit 



sine [ciji^CkRs] . 



In Eq. (U), Hq describes the nonmagnetic (NM) scatter- 
ing region, 



(j^'j In the subsection, we examine the pumping current 

at low frequency limit while maintaining pumping am- 
plitude finite. In this limit, the system is nearly in 
equilibriiUm— and we will use the equilibrium Green's 
functionE^^Ell to characterize the pumping process. Using 
the distribution function, the total charge in the system 
(4) during the pumping is given by 



Vp is the time-dependent pumping potential and Ht de- 
scribes the coupling between electrodes and the NM scat- 
tering region with hopping matrix Tkan- To simplify the 
analysis, we assume the hopping matrix to be indepen- 
dent of spin index, hence 



+ c.c. 



(5) 



In these expressions Cka — ^1 + with a — L,R; c\^^ 
(with a —\, I or ±1 and a = —a) is the creation opera- 
tor of electrons with spin index a inside the a-electrode. 
Similarly d\^^ is the creation operator of electrons with 
spin cr at energy level n for the NM scattering region. In 
writing down Eqs.(|^) and (|^), we have made a simplifi- 
cation that the value of molecular field M is the same fat 
the two electrodes, thus the spin- valve effect is obtainedEiJ 
by varying the angle 0. Essentially, M mimics the dif- 
ference of density of states (DOS) between spin up and 
down electronsEil in the electrodes. In this paper, we only 
consider the single electron behavior. The charge quan- 
tization is not considered so that our system is not in the 
Coulomb blockade regime. In addition, for the nonmag- 
netic regions we are interested in, the Kondo effect can 
be neglected. 

To proceed, |-we first apply the following Bogliubov 
transformatioiiESl to diagonalize the Hamiltonian of the 
right electrodeEj, 



CkRo 



cos{9/2)CkRa - (Tsai{9/2)CkRs 



4r. - cos(0/2)C+,^ - fTsin(0/2)C+,, 
from which wc obtain the effective Hamiltonian 



(6) 



(7) 



kRa 



- cr sin - C|^- ) + c.c] . (9) 

In the following subsections, we will consider two cases: 
(1). parametric pumping in the low frequency limit with 
finite pumping amplitude. (2). pumping in the weak 
pumping limit with finite pumping frequency. 



Q{x,t) - -^q / {dE/2n){G<{E,{V{t)})), 



(10) 



where G< is the lesser Green's function in real space, 
X labels the position, and {V^(t)} describes a set of ex- 
ternal parameters which facilitates the pumping process. 
Within Hartree approximation, is related to the re- 
tarded and advanced Green's functions C and G°, 

G<{E,{V}) - G^E,{V})zJ2t^ME)G%E,{V}) . 

a. 

(11) 

where fa{E) = f(E — qVa). In the low frequency limit, 
the retarded Green's function in real space is given by 



G^{E,{X}) = 



1 



E^Hn 



(12) 



where S'' = J^a is the self energy, and Fq = 
— 2/m[5]^] is the linewidth function. In above equations, 
Q'-.a.< denotes a 2 x 2 matrix with matrix elements G''''lf 

cr,cj 

and a =t, '] . Vp = Vpl where I is a 2 x 2 unit matrix 
in the spin space. In real space representation, is a 
diagonal matrix describing the variation of the poten- 
tial landscape due to the external pumping parameter 
V. The self-energies are givero 



and 



(13) 



(14) 



(8) with the rotational matrix for electrode a defined as 



R 



cos 9^/2 sin 0^/2 
— sin6'a/2 cos 6*^/2 



(15) 



Here angle 9a is defined as 6*^ = and 9r = 9 and EJJ 
is given by 
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and Tl^ 
9^0. 



-2/m(S^g.) is the linewidth function when 
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In order for a parametric electron pump to function 
at low frequency, we need simultaneous variation of two 
or more system parameters controlled by gate voltages: 
Vi{t) — Vio + Vip cos{ujt + . Hence, in our case, the po- 
tential due to the gates can be written as = ^ ■ ViAi, 
where is potential profile due to each pumping poten- 
tial. For simplicity we assume a constant gate potential 
such as that (A.i)x:x is one for x in the first gate region 
and zero otherwise. If the time variation of these param- 
eters are slow, i.e. for V{t) = Vq -I- 6V cos{ujt), then the 
charge of the system coming from all contacts due to the 
infinitesimal change of the system parameter {SV —> 0) 
is 



dQp{t) =J2dv.Tr[Q{x,t)] SV.it) 



(17) 



it is easily seen that the total charge in the system in a 
period is zero which is required for the charge conserva- 
tion. To calculate the pumped current, we have to find 
the charge dQpa passing through contact a due to the 
change of the system parameters. Using the Dyson equa- 
tion dviG^ = CAiC, the above equation becomes, 

/rl W 
— ^TriC-A.-G^r^C^ 

+G'T^G'^A,G'']/0(i?)5F,(t) 

= -1 [ ^EE/^Tr[9s[GT^G'^A,]]<5F,(i) 

where the wideband limit has been taken. Integrating by 
part, we obtain 

dQpit)^q I dEY,idEMY.^5V,it) (18) 
where we have used the injectivityS 



'""^ . ^Tr[G^r,G"A,] 



(19) 



Using the partial density of states dNa0 /dVj defined asB 



^ = -lTr[G^r„G'-A, + c.c.]S^p 
dVj An 

+— Tr[iG'T/3G''r„G''Aj + c.c] 
47r 



(20) 



with dNap/dVj — dNp/dVj, we obtain 
dQp.it) = / dEj2i-dEfp)J2^6V,{t) (21) 

If we include the charge passing through contact a due 
the external bias, then 



dQ.it) = -qf dEY,i~dEf(3)Y. 

-q J di;ETr[r„G'T/3GT(/„-/0)dt 



dNafs dVj it) 



dVj dt 



dt 



(22) 



Furthermore, the total current flowing through contact a 
due to both the variation of parameters Vj and external 
bias, in one period, is given by 



dt dQa/dt 



(23) 



where r = 27r/a; is the period of cyclic variation. If there 
are two pumping parameters, Eq.(]2^) can be written as 
when a = 1, 

T J J . 

dN II dV jit) dN 12 dVj it) 



dV-i dt 



dV, dt 



]dt 



-^JdtJ di;Tr[riG'T2GT(/i -/2) 



(24) 



and 



2t 



dt 



dNii dVjit) 



dE dEih 
dNi2 dVjit). 



dV, dt 



dt 



where Ji = + J)^' ■ In Ref. JY' has been iden- 
tified as the current due to the external bias and J^^-* as 
pumping current. For = 0,7r, all the 2x2 matrix are 
diagonal. In this case, it is easy to show that Eq.(|2^) 
agrees with the result of Ref. |2^. If the external bias is 
zero, Eq.(^) reduces to the familiar formulati when there 
are two pumping potentials. 



r(2) 



dt 



(1) 



(25) 



quj 
2^ 



dt 



dNa dXi dNa dX2 



dXi dt dX2 dt 



(26) 



B. Finite frequency pumping in the weak pumping 
limit 



In this subsection, we will calculate the pumping cur- 
rent at finite frequency. The Keldysh nonequilibrium 
Green's function approach p^ed here is in the standard 
tight-binding representationEj. We could not use the mo- 
mentum space version because the time dependent per- 
turbation (pumping potential) inside the scattering re- 
gion is position dependent. Hence it is most suitable 
to use a tight-binding reaLspace technique. In contrast, 
in previous investigationsEa on photon-assisted processes 
the time-dependent potential is uniform throughout the 
dot and therefore a momentum space method is easier to 
apply. 

Assuming the time-dependent perturbations located at 
the different sites, j = io,jo, and fco etc, with 



Vjit) = Vj cosiujt + 4>j) 



(27) 
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When there is no interaction between electrons in the 
ideal leads L and R, the standard nonequilibrium Green's 
function theory gi^fes the following expression for the time 
dependent currenttJ, 

Jc.{t) = -q f diiTr[G'-(i,ti)S<(ti,i) 

J — OO 

+ G<(t,ti)S^(ii,t) + c.c.] (28) 

and transmission coefhcient 

T{E) = Tr[ri{E))G'-{E)TUE))G^{E)] (29) 

where the scattering Green's functions and self-energy 
are defined in the usual manner: 

G:';..'ih,h) = T^^{±tlTt2){{dUtl),d+^,{h)}) (30) 
G< .,..(^1,^2) = i{d+,,{t2)d,..Xti)) (31) 



(32) 



The average current JL{t) from the left lead can be 
written as 

<JL{t)>^--f dt f dtiTr[Gl,{t,h)i:<{h,t) 

+ G<(i,ti)S2(ti,t)+c.c.] (33) 

In the absence of pumping, the retarded Green's func- 
tion is defined in terms of the Hamiltonian iJo, 



G"''(£;) 



1 



(34) 



and G"^ is related to the retarded and advanced Green's 
functions G^'' and G"'' by, 



G"<(i;) = G"''(-B)S<(i;)G"''(£;) 



(35) 



Now we make use of the time-dependent pumping po- 
tentials Vj{t) as the perturbations to calculate all kinds 
of Green's functions up to the second order, and corre- 
sponding average current. 

First, we calculate the current corresponding to the 
term G\i{t^ti)Yif^{ti,t) in Eg. (p3|) . Dyson equation for 
G\i{t,ti) gives the second order contribution: 

G?i\tM) [ [ d^dyG%{t - x) 



V,{x)G%{x-y)V,{y)Gf,{y^h) 
^^G°,^V,G°^VfeGri 



(36) 



Substituting Eq.(|36|) into Eq.(p3|) and completing the 
integration over time x, y, ti, t, it is not difficult to calcu- 
late the average current < Jli > due to the first term in 
Eq.(^ (see appendix for details). 



jk 

[G%iE.)e^'^''^ + G%iE+)e'^^''^]Gll] (37) 



where Akj = 4>k — (t>j is the phase difference, E± ~ E±uj, 
and G''''''< = G''''''<{E). 

Now we calculate the second term Gf;^(t, ii)S|^(ii, t) 
in Eq.(ll). Using Keldysh equation, G< = G''S<G'^ , 
we have 



G< f-^r ■sr^<f~ia , f-^r \^<{^a 
jk — ^il^L^lk + ^jN-^R^A 



Nk 



(38) 



where S< = iTafa- Expanding G^'"' up to the second 
order in pumping parameters, we obtain the second order 



contribution from G^-^, 



G{2)< f~i7' ^<f~ia I i^r \y</^a 
11 — *^11^L*-^11 + '^ITV^fl'^Ari 

— "^11 -^L^^ll+^ll ^L*^ll +^ll^L*Jll 

+ ^1N ^R^Nl + ^IN ^R^Nl +^1N^R^N1 



= ^[G?^^V,G0^V,.G0.< + G;,^V,G°<V,G0^ 



Jk 



G?<V,GO,"VfeG°?] 



(39) 



where we have used Eq. (|38[) to simplify the expression. 
After some algebra, we have the following three expres- 
sions corresponding to each term in Eq.(|39|), 



E 

jk 



qVjVk fdE 
2tt 



-Tr S? G 



G%iE+)]G 



kl J 



(40) 



ok 



\Gti 



(41) 



qV.Vk 



jk 



r'0< 
l'^ij 



[G%iE^)e^^''^ 



(42) 



The final pumped current is the sum of Eq. (|37|) , ( |4C. , 
(^l|), ( ^2|) and their complex conjugates, in addition to 
the current directly due to the external bias (see second 
term of Eq.(|^). If the external bias is zero, the expres- 
sion of the pumped current can be simplified significantly. 
Note that in the equilibrium, the lesser Green's function 
satisfys the fluctuation-dissipation theorem. 



G°<{E) = -f[E) [G0'^(£;) - G'^'^iE)] 



(43) 



and 



G"<(i?±) = -f{E±) [G'"^(i?±) - G^^E^)] (44) 
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(11 + ®+®* leads to 
qVjVk 



jk 



J f Tr [r,G?;:/(i?) 



while @* + ® *+(||) gives to 



(45) 



lOa/ 
'jk 



[G°liE.)e^^''^ + e~^^''^G°UE+)]Glt] 



(46) 



furthermore, (^T])+c.c. becomes 

^ 7,. 



[/_(G,"^(i?_)-G"-(i?_)) 



/+(G°^(ii;+)-G^",«(i?+))]G°?] 



(47) 



where /± = f{E±). 

Combining Eqs.(^), (^) and (^), we finaUy obtain 



qVjVk 



ife 

[(/_-/)(GO;:(i?_)-G!^,"(i?_))e^^'=^ 

+e-^^^(/+ ^ /)(G0^(£;+) - GO,"(i?+))]GO^] (48) 

In the hmit of small frequency, we expand Eq.(^) up to 
the first order in frequency and use the fact that 



and Dyson equation, 



GOr/-iOr 
ji ^ik 



dG% 



(49) 



(50) 



we obtain, 



qujVjVk sm{Akj) f dE 



2n 



dsfiE) 



Tr{r,.G;^-(ii;) [G%iE) ~ Gf.iE)] GlUE)} 



Jk 



Tr 



dG^ aG?? 



J 27r 



(51) 



which is the same as Ref. 0whcn 9^0. 



III. RESULTS 



We now apply our formula Eqs. 



a TMR junction. For current perpendicular to the p 



1) and (H) to 



ane 



geometry, the TMR junction can be modeled by an one- 
dimensional quantum structure with a double barrier po- 
tential U{x) = Xi6{x + a) + X2S{x — a) where 2a is the 
well width. For this systeipjbhc Green's function G{x, x') 
can be calculated exactl5/c3. The adiabatic pump that 
we consider is operated by changing barrier heights adi- 
abatically and periodically: Xi — + Vp sin(cL;t) and 
X2 = Vq -\-Vp sin(wt + (f). In the following, we will study 
zero temperature behavior of the pumped current. In 
the calculation, we have chosen M = S2.0 and Vb = 79.2. 
Finally the unit is set by ?i = 2m = Icll. 

We first study the pumped current with two pump- 
ing potentials in the adiabatic regime. Fig.l depicts the 
transmission coefficient T versus Fermi energy at several 
angles 9. As expected, among different angles, T is the 
largest at = and smallest a± 9 — it with the ratio 
Tmax^ lTmax{T^) ^ 4. This givcs the usual spin valve 
effectEil. Fig. 2 plots the pumped current versus Fermi 
energy at different 9. Here we have set the phase dif- 
ference of two pumping potentials to be 7r/2. Similar to 
the transmission coefficient, we obtain largest pumped 
current at = and smallest current at = tt. We 
found that the ratio ImaxiO) / Imax(j^) is about the same 
as that of transmission coefficient. As the pumping am- 
plitude doubles, the peak of pumped current is broadened 
and the maximum pumped current is doubled (see inset 
of Fig. 2). This is understandable since at large pump- 
ing amplitude the instantaneous resonant level oscillates 
with a large amplitude and hence can generate heat cur- 
rent in a broad range of energy. The spin valve effect of 
pumped current is illustrated in Fig. 3 where the pumped 
current versus 9 is shown when the system is at reso- 
nance. We see that the pumped current is maximum at 
9 = and decreases quickly as one increases 9 from to 
TT. For larger pumping amplitude, we have similar behav- 
ior (inset of Fig. 3). In Fig. 4, we plot the pumped current 
as a function of phase difference (fi between two pumping 
potentials. We see that the pumped current is antisym- 
metric about the (f> = n. The nonlinear behavior is clearly 
seen which deviates from the sinusoidal behavior at small 
pumping amplitude. In Fig. 5, we show the pumped cur- 
rent in the presence of external bias. In the calculation 
we assume that Vl = —ojV/2 and Vr = ujV/2 so that 
the external bias is against the pumped current when V 
is positive. Due to the external bias, the total pumped 
current (dashed line) decreases near resonant energy and 
reverses the direction at other energies. Now we turn to 
the case of finite frequency pumping. We first present 
our results (Fig. 6 to Fig. 8) at small pumping frequency 
uj = 0.002. In Fig. 6 we plot the pumped current as a 
function of phase difference (f) near the resonant energy. 
At = 0, the magnitude of pumped current is much 
larger than that at 6' = 7r/2 or tt. We notice that at 
= and (f) = TT, the pumped current is nonzero similar 
to the experimental ajiomaly observed experimentally for 
nonmagnetic systemH. The pumped current away from 
resonant energy is shown in Fig. 7. At = 0, we see 
that the pumped current is sharply peaked at resonant 
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energy. The pumped current is positive for = 7r/2 and 
negative for = and tt. At 9 = 7r/2 or tt, the pumped 
current at = 7r/2 (dashed Une) is much larger than that 
at other angles. Fig. 8 displays the pumped current as a 
function of 9 near resonant energy. For cj) = tt/2 (dotted 
line), we see the usual behavior that large pumped cur- 
rent occurs at 9 = and it decreases to the minimum at 
9 = IT. For = or TT, however, we see completely differ- 
ent behavior. The pumped current is the still the largest 
at 9 = but the direction of the pumped current is re- 
versed. As one increases 9, the pumped current decreases 
and reaches a flat region with almost zero pumped cur- 
rent. Now we study the effect of frequency to the pumped 
current versus 9 (see Fig. 9). We will fix the phase differ- 
ence to be = 7r/2 and energy near resonance. At small 
frequency lu = 0.002 (dashed line), the pumped current 
versus 9 show usual behavior. When the frequency is in- 
creased to lli = 0.004 (dot-dashed line), two peaks show 
up symmetrically near 9 = ±7r/4 while the minimum 
is still at 6* = TT. As frequency is increased further to 
u) = 0.006 (dotted line), the pumped current near 9 = 
reverses the direction and the new peak position shifts 
to ~ O.SStt. Upon further increasing w, the curve of 
pumped current versus 9 develops a flat region between 
9 = ±0.357r with positive current while the magnitude 
of the negative pumped current at 9 = becomes larger 
(see Fig. 9b). Finally, at even larger frequency uj = 0.1, 
all the pumped currents are negative. This behajdpr can 
be understood from the photon assisted procesaij. The 
quantum interference between contributions due to pho- 
ton emission (or absorption) near two pumping potentials 
is essential to understand the nature of pumped current. 
In addition to the interference effect, the pumped current 
is also affected by a competition of between the photon 
emission and absorption processes which tend to cancel 
to each other. It is the interplay between this compe- 
tition and interference that gives rise to the interesting 
spin valve effect for the pumped current. Finally, we show 
in Fig. 10 the pumped current versus pumping frequency. 
We see that at small frequency the current is positive and 
small, at large frequency the current is much large and 
is negative. 

In summary, we have developed a general theory for 
parametric pumping in the presence of ferromagnetic 
leads. Our theory is based on the nonequilibrium Green's 
function method and is valid for multi-modes (in two or 
three dimensions) and can be easily extended to the case 
of multi-probes (although most of calculations are for 
two probes). In the parametric pumping, two kinds of 
driving forces are present: multiple pumping potentials 
inside the scattering system as well as the external bias in 
the multi-probes. Two cases are considered. In the adi- 
abatic regime, the system is in near equilibrium. In this 
case our theory is for general pumping amplitude. At fi- 
nite frequency, the system is away from equilibrium. Our 
theory is up the quadratic order in pumping amplitude. 
This theory allows us to examine the pumped current in 
broader parameter space including pumping amplitude. 



pumping frequency, phase difference between two pump- 
ing potentials, the angle between magnetization of two 
leads. 
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IV. APPENDIX 



Now we show that 



B = -f dt [ dtiTr[Fo(ii,t)G(t,ti 



where 



'F2{E_)]F3iE)] 



G = FiV.FaV^Fa 



(52) 



(53) 



and F,'s {i = 0,1,2,3) satisfy F,{ti,t2) = F,{ti - is). 
Taking the Fourier transform. 

Fit) = ^ J dEe-'^'F{E) 



(54) 



we obtain 



B=iJ\tj' dhj ^Fo(i;)e-'^(*^-*) J dxdy 
[Flit - a;)F4a;)F2(x - y)Vpiy)F3iy - t^)] 

^ nf/S/:*/-. 

FoiE)FiiEi)V^iE2)F2iE3)Vf!iEi)F3iE5) 

^i(E-Ei)t^i{E^-E)t^^i{Ei-E2-E3)x^i{E3-Ei-E^)y 

where 

Yo^iE) = T:Vc,[e'^-5iE+) + e-**°^(£;_)] (55) 
Integrating over x and y yields, 

I r , f -n dE, f dE f* , 

Fi(Si)V„(£;2)F2(^3)V^i(S4)F3(S5)e'^^~^^^* 
g»(iJ5-s)t,(2;r)25(£;i -E,- E;)SiE, - E^ - E,) 

Integrating over ti,t and using Eq.(|5a), we have 
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B = 



dE2dEidE5 f dE Fo{E) 



2-K 27r 2-K J 2-K i[E^ - E - i5] 
F^iE^ +Ei + Ez)\^{E2)F2{Ei + E^) 
Yp{Ei)¥^{E^)5{E2+E^) 
_ Vc^Vp fdE,dE FojE) 

4 J 271 2tt i[E5 - E ~ iS] ^'^ 

[F^iEr, + w)e^^'^° + e~''^^-F2{E5 - u;)]F3iE5) 
Using the theorem of residue, we have 

'^^ ■iw^'^f ■. =^^Po{E5) (56) 
i[E5 — E — id\ 

we thus obtain Eq.d53). 
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FIG. 1. The transmission coefHcient as a function of Fermi 
energy at different angle 6 between the magnetizations of two 
leads: 9 = (dashed line), 6 = -k/2 (solid line), and 6 = 1: 
(dotted line). 

FIG. 2. The pumped current as a function of Fermi energy 
at different 0: 6 — (solid line), 9 — n/2 (dot-dashed line), 
and 9 = -IT (dotted line). Other parameters are (/> = 7r/2 
and Vp = O.OSVb- Inset: the same as the main figure except 
Vj, = 0.1Vo. 

FIG. 3. The pumped current as a function of 9. Here 
4> = 7r/2, Ef = 37.55 and Vp = 0.05Vb. Inset: the same as 
main figure except Vp = O.lVb. 

FIG. 4. The pumped current as a function of phase dif- 
ference (f> at different 9. 9 = (dashed line), 9 = it/2 
(dotted lino), and 9 = n (solid line). Other parameters are 
Ef = 37.55 and Vp = 0.05Vo. 

FIG. 5. The pumped current as a function of Fermi energy 
with an external bias Vl — Vr = —0.02a; at different 9: (a). 
9 = IT. (b). 9 = n/2. (c). 9 = 0. In (a), (b) and (c), solid line: 
pumped current, dotted line: current due to external bias, 
dashed line: total current. Here = tt/2 and Vp = 0.05Vo. 

FIG. 6. The pumped current as a function of 4> at finite 
frequency. Dashed line: ^ = 0, dotted line: 9 = it/2, solid 
line: 9 = it. Here Ef = 37.55 and lj = 0.002. 



FIG. 7. The pumped current as a function of Fermi energy 
at finite frequency at different 9: (a). 9 = 0. (b). 9 = tt/2. 

(c). 9 — n. In (a), (b) and (c), dotted line: (j> = 0, dashed 
line: 4> = n/2, solid line: (j) = n. Here w = 0.002. 

FIG. 8. The pumped current as a function of 9 at finite 

frequency. Here dashed lino: (j) = 0, dotted lino: 4> = n/2, 
solid line: 4> = n. Other parameters are Ef = 37.55 and 
uj = 0.002. 

FIG. 9. The pumped current as a function of 9 at dif- 
ferent frequencies. (a), cu = 0.002 (short dashed line), 
0.004 (dot-dashed line), 0.006 (dotted line), 0.008 (solid line), 
(b). UJ = 0.01 (solid line), 0.02 (dot-dashed line), 0.05 
(short dashed line), 0.1 (dotted line). Other parameters are 
Ef = 37.55 and 4) = tt/2. 



FIG. 10. The pumped current as a function of frequency. 
Here 9 = 0,(f> = n/2 and Ef = 37.55. 
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